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We analyze the possibility of delensing Cosmic Microwave Background (CMB) polarization maps 
using foreground weak lensing (WL) information. We build an estimator of the CMB lensing poten- 
tial out of optimally combined projected potential estimators to different source redshift bins. Our 
estimator is most sensitive to the redshift depth of the WL survey, less so to the shape noise level. 
Estimators built using galaxy surveys like LSST and SNAP recover up to 80-90% of the potential 
fluctuations power at I < fOO but only «10-20% of the small-angular-scale power (I < 1000). This 
translates into a 30-50% reduction in the lensing B-mode power. 

We illustrate the potential advantages of a 21-cm survey by considering a fiducial WL survey 
for which we take the redshift depth z max and the effective angular concentration of sources n as 
free parameters. For a noise level of 1 fiK arcmin in the polarization map itself, as projected 
for a CMBPol experiment, and a beam with #fwhm=10 arcmin, we find that going to z max =20 
at n — 100 galaxies/arcmin 2 yields a delensing performance similar to that of a quadratic lensing 
potential estimator applied to small-scale CMB maps: the lensing B-mode contamination is reduced 
by almost an order of magnitude. In this case, there is also a reduction by a factor of ~4 in the 
detectability threshold of the tensor _B-mode power. At this CMB noise level, the B-mode detection 
threshold is only 3x lower even for perfect delensing, so there is little gain from sources with 
Zmax > 20. The delensing gains are lost if the CMB beam exceeds ~ 20 arcmin. The delensing gains 
and useful z max depend acutely on the CMB map noise level, but beam sizes below 10 arcmin do not 
help. Delensing via foreground sources does not require arcminute-resolution CMB observations, a 
substantial practical advantage over the use of CMB observables for delensing. 



PACS numbers: 98.62.Sb, 98.80.-k, 98.70.Vc, 98.80.Es 



I. INTRODUCTION 



The anisotropies in the CMB have been long recog- 
nized as a major probe for cosmology. The WMAP 
satellite has measured the temperature fluctuations up 
to l m ax <=1000 and has confirmed the cosmological 
standard model of a power-law, flat, ACDM universe. 
Much hope lies with CMB polarization measurements 
because they have the potential to unveil some of the 
unknowns of inflation. While some predictions of infla- 
tion (such as nearly flat space curvature, nearly scale- 
invariant power spectrum, and Gaussianity of the primor- 
dial fluctuations) have been confirmed by CMB and large 
scale structure experiments, there is another prediction 
which is yet to be verified. This is the existence of an 
almost scale-invariant spectrum of gravitational waves, 
whose amplitude is directly related to the energy scale of 
inflation. The inflationary gravitational waves are ten- 
sor perturbations to the metric and we expect them to 
leave a curl-like signature in the CMB polarization field, 
i.e. a B-mode pattern. Density fluctuations, arising from 
scalar perturbations to the metric, create a gradient-like 
component in the polarization field, the E mode. In 
the linear regime, density fluctuations do not create a 
B mode, so a detection of the latter would confirm the 
existence of primordial gravitational waves; it would also 
provide the energy scale of inflation, which we could use 
to distinguish between different inflationary scenarios. 

CMB polarization measurements are difficult to carry 
out, primarily for two reasons. Firstly, the amplitude 



of the signal is very small: for example, the scalar E- 
mode power is 2-3 orders of magnitude smaller than the 
scalar temperature power, tensor i?-mode much lower. 
Secondly, the polarization foregrounds are very poorly 
understood and they dominate the CMB signal at al- 
most all relevant frequencies. Therefore, foreground re- 
moval and detector sensitivity are two of the most stifling 
limitations of a polarization experiment. 

£>-mode measurements are additionally obstructed by 
WL contamination, especially of the recombination sig- 
nal (1<100). Thus CMB polarization measurements are 
able to provide inflationary insights to the extent to 
which we can: 1. remove the foreground contribution to 
the overall signal and 2. delens the B mode power and 
extract the tensor contribution to it. The delensing pro- 
cess consists of the reconstruction of the lensing potential 
from chosen observables. The estimated lensing potential 
is then used to evaluate the WL-created -B-mode signal 
and to subtract it from the measured i?-mode map. 

In this paper we probe the ability of weak lensing (WL) 
galaxy surveys to delens the CMB in the absence of fore- 
grounds. To be specific, we try to answer two ques- 
tions: what attributes should a galaxy or 21 cm WL 
survey and also a CMB polarization mission have in or- 
der to detect the tensor B mode? What is the mini- 
mum amplitude of the B mode, expressed in terms of 
the tensor-to-scalar ratio, r, that can be detected using 
galaxy or recombination observations for delensing? We 
take three examples of surveys to illustrate how our es- 
timator works: the ground-based Large Synoptic Survey 
Telescope (LSST |2g|), the space-based Supernova Accel- 
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eration Probe (SNAP [231) and a toy model mimicking 
recombination-era 21-cm observations that we mention in 
more detail in section mv\ The outline of the paper is as 
follows. In SjTT]we describe the WL contamination of the 
tensor B mode. In i]IIII we present our minimum-variance 
lensing estimator. In EHVI we determine the minimum de- 
tectable r when we use this estimator to delens. In <|V] 
we discuss our results and draw conclusions. Let us now 
briefly mention similar work existing in the literature. 
There has been a vast and impressive amount of work 
on the topic of CMB delensing and r-detection. The 
great majority of this work uses the CMB observables 
6, E, B to reconstruct the projected potential. Averag- 
ing over various quadratic combinations of the tempera- 
ture field (e.g. see the work of Zaldarriaga and Seljak [l|, 
Bernardeau [2j , Hu [H, Hu [1]) and the polarization field 
(e.g. Guzik et al. @, Hu and Okamoto Q, Kesden et al. 
p}) has been thoroughly considered for the projected po- 
tential reconstruction. To give a quick summary: one can 
build minimum- variance, unbiased estimators, using cer- 
tain field statistics, as shown by Hu [3] ■ For a post-Planck 
experiment (sensitivity of 0.3 /^K arcmin and beam size 
of 3 arcmin), Hu and Okamoto tj found that the most 
efficient of these estimators can map the potential up to 
K1000. Kesden et al. § and Knox and Song used this 
last estimator to predict the minimum detectable r as a 
function of CMB experimental characteristics. There is 
another, more promising method for lensing potential re- 
construction, based on likelihood techniques. Hirata and 
Seljak jlOj have built a maximum-likelihood estimator for 
the convergence field using temperature maps and have 
found its performance similar to that of the quadratic 
estimator introduced by Hu 0]. Hirata and Seljak [llj 
found that the same maximum likelihood estimator built 
from polarization maps is even more effective: there is an 
order of magnitude reduction in the mean squared error 
in the lensing reconstruction compared to the quadratic 
estimator method, if the survey characteristics are ade- 
quate: sensitivity of 0.25 fiK arcmin and a beam size of 
2 arcmin. Amarie et al. 12| analyze the detectability of 
tensor B modes in the presence of polarized dust emis- 
sion, as a function of sky coverage; Verde et al. [l3[ do a 
study of optimal surveys for i?-mode detection, consid- 
ering both dust and synchrotron emissions. One disad- 
vantage of the reconstruction methods presented so far 
is that they require high-resolution CMB maps: they use 
arc-minute structures of the CMB fields, to reconstruct 
degree-scale maps of the deflection field, as explained by 
Hu [J] . Sigurdson and Cooray [14[ point out that one can 
use non-CMB observables to delens the CMB; in this case 
the requirement for high angular resolution of the CMB 
mission can be relaxed significantly. These authors de- 
termine lower limits for the detectable r using the 21 
cm radiation emitted by neutral hydrogen atoms to de- 
lens the CMB. We follow a similar approach here, but 
employ foreground galaxies instead of 21-cm emission as 
the source plane for delensing. 



II. B MODES AND WEAK LENSING 

As mentioned in fJH inflation predicts scalar and ten- 
sor perturbations; the former are responsible for density 
fluctuations and for the formation of structure. The lat- 
ter are gravitational waves. Direct detection of gravita- 
tional waves has been a physicists' dream for a long time 
and there is ongoing effort to make this dream come true 
(e.g. see LISA |2j|, LIGO [IH). There is however another 
means to ascertain the existence of tensor perturbations, 
by measuring the CMB polarization. 

The E and B modes are defined as linear combina- 
tions of the Stokes parameters Q and U, such that the 
power of B arising from linear order scalar perturbations 
is zero, as shown by Zaldarriaga and Seljak [l5[. Tensor 
perturbations yield both E and B modes, roughly of the 
same magnitude; both modes could in principle provide 
information on gravitational waves and inflation, but in 
practice the tensor E modes are overwhelmed by their 
scalar counterparts. The amplitude of tensor modes is 
quantified by the tensor-to-scalar ratio r, the ratio of the 
expectation values for the quadrupole tensor and scalar 
temperature anisotropies. 

Even if one assumes flawless detector technology, B- 
mode measurements are plagued by the presence of po- 
larized foregrounds and by WL contamination. A de- 
tailed discussion of the foreground problem is beyond the 
scope of this work and we refer the interested reader to 
the Final Report of the Task Force on CMB Research 
[ill. We briefly mention that while most galactic polar- 
ized foreground emission is currently not well known at 
the frequencies of interest to CMB experiments, there is 
reason to believe that this situation will change in the 
near future. Indeed, the WMAP satellite is gathering 
all-sky data on the synchrotron emission and the HFI in- 
strument on the Planck satellite will map the dust emis- 
sion. Numerous smaller, ground-based missions, with 
higher angular resolution, plan to complement the above- 
mentioned satellite data. All this information should en- 
able the study of polarized foreground emission as a func- 
tion of both frequency and angular scale. 

We now sketch the equations that describe the WL 
contamination of the tensor B mode. The groundwork 
on this subject was done by Zaldarriaga and Seljak 171: 
for a more recent review see Lewis and Challinor |18 |. 



A photon travelling from the last scattering surface is 
affected by potential wells along the line of sight. Its 
frequency changes because the potentials vary with time 
(the Rees-Sciama effect, the integrated Sachs- Wolfe ef- 
fect) and its path can be transversally shifted by lensing. 
To a good approximation, the lensed wave has the same 
Stokes parameters along a path n as the unlensed one, 
only shifted by the average deflection along that path: 




where we denote by X the lensed field of parameter X, 



3 



and a is the deflection angle. To linear order, the de- 
flection angle is the gradient of the projected potential: 
a=Vip. Assuming that deflection angles are very small 
compared to the characteristic scale of fluctuations, we 
can Taylor expand the lensed parameters around the un- 
deflected path x. Using the flat-sky approximation, we 
can write the lensed B mode in Fourier space as: 



d 2 V 



B{l)=B{l)- I -^rf(l',l) HL)E(l') 



d 2 V d 2 V 
~2tt~~2tt 



■g(l',l",l)^(l")r(l"-L)E(l'). (2.1) 



Here B denotes the tensor B mode and E is the unlensed, 
scalar E mode. The above equation neglects the lensing 
effect on the tensor E and B fields, since it is very small 
compared to the lensing of the scalar E mode. The func- 
tions / and g are defined by: 

/(I',Z) = r-Zsin[2((fc/-^)] 
g {V, I", l)=l>- (I"- L){V ■ I") sin[2(^ - ft)] , (2.2) 

where L = 1 — 1' and (pi is the angle between I and the x 
axis. We assume statistical isotropy and parity invariancc 
and we adopt the power spectrum definition of Lewis and 
Challinor O: 



(E(l)E*(l')) = 5(l-l>)CF, 

J l > 



(B(l)B*(l')) = 6(l-l')C» 



(B(l)E*(l)) = 0. 



(2.3) 



With the above considerations in mind, we write down 
the lensed B-mode power spectrum as a sum of the ten- 
sor (primordial) B-mode power and the weak lensing B- 
mode power: 



with the appropriate definitions: 



£iB ^jtensor 

cr = J £L\fV,i)\ 2 ctc$+ 



(2.4) 



(2.5) 



(27T) 

i r d 2 V d 2 i" 

'I J (2^)2 (2tt) 2 



\g{V,l",l)?C%Cl„_ L ^, (2.6) 



where C Z B is the tensor B mode power spectrum, Cf is 
the power spectrum of the unlensed scalar E mode and 
Cf is the lensing potential power spectrum. 



space of the scalar E-mode power with the lensing poten- 
tial power, as shown in Eq. (|2.4[) . In order to have access 
to the tensor B mode, we need to subtract the lensing 
contribution from the measured B-mode map. The im- 
pediments are that we know neither the lensing potential 
nor the unlensed E mode power. 

In this section we address the first of these prob- 
lems and build a lensing potential estimator using WL 
foreground galaxy survey information, independent of a 
CMB experiment. Unlike the estimators built out of 
CMB observables, our estimator does not require ex- 
tremely high resolution in the CMB polarization map. 
This is an important advantage, since high resolution is 
the main cost driver of CMB polarization experiments. 
The disadvantage of our approach is that WL surveys de- 
tect source galaxies up to a certain redshift, depending on 
the specifics of the survey. Our estimator cannot recon- 
struct the potential fluctuations evolving between that 
upper-limit redshift and the redshift of recombination. 

We build the lensing potential estimator as a weighted 
sum of the projected potential measured to different 
source galaxy redshift bins. We choose the weights so 
that the estimator be optimal, i.e. the variance of the er- 
ror in each mode is minimal. We work in the flat sky ap- 
proximation and we assume a flat Friedmann-Robertson- 
Walker (FRW) metric. We write the estimator as: 



(3.1) 



where the sum is over source bins i and the a's weight the 
contribution of the projected potential to each redshift 
source bin. The projected potential to the ith source bin 
is defined by: 



MO) 



dxWi(xW u (xe,x), (3-2) 



where $ 3Z? (x<?,x) is the 3D gravitational potential spec- 
trum. In Eq. (|3.ip we use the Fourier transform of the 



projected potential: ipi(l) = J d 2 



il ■ e 



MO)- The 



lensing weight function Wi (x) °f the- source bin i is given 
by: 

r;; +1 dz s v(z s ) x ty- x u( Zs - z{x)) 

lW c 2 J z Zi+1 dz s P(z s ) 

(3.3) 

Here u(x) is the Heaviside unit step function and V(z s ) 
is the redshift distribution of the source galaxies. The 
residual of our estimator, defined as: 



III. THE CMB POTENTIAL ESTIMATOR 



A. The ideal case 



tells us how the estimator compares to the true CMB 
projected potential: 



The measurable .B-mode power is a sum of the primor- 
dial B-mode tensor power and a convolution in Fourier 



i> C un(0)= d x w CMB (x)$ 3D (xO,x), 

Jo 
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where W CMB (x) is the same weight as in eq. (|3.3[) . but 
written for the particular case of one source at rcdshift 
z C mb , i-e. V{z s ) = S D (z CMB - z s ): 

, s _ 2 x(z CMB )-y 
W CM n[X) ~ "a / s • 

C Xl Z CMB,iX 

We determine the weights by requiring that the vari- 
ance of the above-mentioned residual be minimal: 



of galaxies to limit the performance of our estimator 
quite significantly. The intrinsic ellipticities of galaxies 
increase the variance of the projected potential in each 
source bin. The shape noise effect on W SS (Z) is: 



d 



da k (l) 



(I^ombCO-^cmbCOI 2 ) =0, Vfc. 



(3.4) 



To see the solution of this equation, let us express the 
variance of the estimator's residual as a matrix equation: 

( |^omb(0 - Vw(0l 2 ) = ■ w ss (z) • a*(t) 

-2a(l) • W SCMB (0 + W CMB CMB (/) . (3.5) 

In the above equation we have used the following defini- 
tions: 

w? OMB (0 = (iM0i£mb(0>. 

W CMBCMB (/) = & CMB (l)^* MB (l)). (3.6) 

Using the relation P$(fc) = [| Hg (1 + z)] 2 X P s {k), 
to go from the potential power spectrum to the density 
power spectrum and also the Limber approximation, we 
calculate the W-matrices: 

Wff(Z) = ^fj^dz^il + zfW^W^Ps^l 

^fr d A 1+z f Wl (x)w CMB (x)Ps(-), 

I Jo dz x 



(0 



^°°d^(l + z) 2 TU c 2 MB ( X )F,(i) 



CMB CMB 



(0 



(3.7) 



where W = 

We are now able to write down the weights a that 
minimize the variance of our estimator's residual: 



a(Z) = W SCMB (0- (W 58 )" (*)• 



(3.8) 



We note that the weights depend only on the magnitude 
of I. Our next goal is to rewrite the expression of the vari- 
ance of 11(1), using eq. (|3.8|) ; the Z-dependence is implicit 
here: 



Var TZ^ = W CMB CMB - W gCMB • (W ss ) _1 • W gc 



(3.9) 



The modes where Var7£(Z) — + are well reconstructed by 
our estimator. 



B. The effect of Shape Noise 

So far we have considered an ideal case, with clean 
WL measurements. In fact we expect the shape noise 



ol 4 



(wff) SJ v(0 = wf?(0 + ^^^, 



(3.10) 



where cr 7 is the uncertainty in the measurement of one 
galaxy and fii is the number density of sources in bin i. 



C. The CMB Lensing Potential Estimator for 
Fiducial Surveys 

We use three examples of WL surveys to probe the ef- 
ficiency of our estimator: (i) the ground-based LSST; (ii) 
the space-based SNAP; (iii)a toy model, which could be 
a future cosmic 21-cm radiation survey. Depending on 
the instrument used, four characteristics describe a WL 
survey: the redshift distribution of the source galaxies, 
"P(z), the angular concentration of source galaxies, n , 
the area of the sky covered, A, and the noise level of the 
intrinsic ellipticities of galaxies, cr 7 . In all cases, we take 
<7 7 =0.3. For observations of 21-cm emission at z > 5, 
the lensing signal might be extracted from shape or den- 
sity correlations of discrete objects (e.g. "minihalos," , 
Barkana and Loeb [23]) while at higher redshifts the de- 
flection field may be derived from an analysis of the full 
3D data cube in a manner analagous to the estimation 
of len sing from the CMB itself (e.g. Zahn and Zaldar- 
riaga 24]). We retain the notation of ct 7 and n even in 
the absence of discrete sources, since the quantity a^/n 
still has meaning as the noise level of the lensing recon- 
struction. The attainable shape noise levels from 21-cm 
observations of z > 5 are highly dependent on the level 
of foreground contamination and the history of reioniza- 
tion, so we will retain the simple n parameterization in 
this paper. For LSST and SNAP we employ the source 
distribution introduced by Smail et al. [l9(: 



(3.11) 



We normalize this distribution so that 
dzV(z) = h . 

The mean redshift of this distribution is (z) = 3zq, where 
zq is a survey parameter. For LSST, we take zq = 0.33, 
n =30 galaxies/ arcmin 2 and .4=20000 deg 2 . For SNAP, 
zq = 0.5, n =100 galaxies/arcmin 2 and .4=1000 deg 2 . 
For the toy model we adopt a "box" (Heaviside) source 
galaxy distribution, i.e. all source bins have equal num- 
ber of source galaxies: 



V(z) =u(z n 



z). 



(3.12) 



Since we are trying to find an optimal WL survey for our 
projected potential estimator, we consider z max and n as 
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free parameters for the toy model: z max ranges from 1 
to 100 and n from 4 to 250 galaxies/arcmin 2 . However, 
we do fix „4=20000 deg 2 , just like for LSST; this seems 
a reasonable choice, which also allows us to compare the 
delensing efficiencies of the two surveys. 

We now specify a few noise definitions for CMB po- 
larization experiments, which will be relevant in section 
E )IVI We adopt the conventions of Knox [2(j. We as- 
sume that different pixels in the polarization map have 
uncorrelated noise with uniform variance o 2 ix . For one 

2 

detector, cr 2 = , where s is the sensitivity of the 

P Ipix 

detector and t P i X is the time spent on the pixel. Then 
the noise power spectrum is: 



C, 



N _ 



N N* 



.)=w- l {W* eam )-' A 5 w 5, 



(3.13) 



Here, w is the weight per solid angle, given by: = 
a pix ®>pix ■ The weight per solid angle has the advantage 
that is independent of the beam size for fixed survey 
time and survey area. We choose the beam to be Gaus- 
sian, therefore the beam window function is: wf am = 
expH 2 a 2 eQm /2), where a 2 beam = 2 fwhm /8ln2. 

We define a weighted estimator for the measured E 
mode, with weights chosen so that the variance of the 
residual of this estimator is minimal. In this case the 
estimator E(l) is just a scalar weight af times the noisy 
observation of E(l). Following the same steps as in the 
case of the lensing potential estimator, the weights and 
the variance of the residual of the measured E mode es- 
timator are: 



of 



Cf 



Cf 



Cf 



V&r1l E (l) 



CfCf 



Cf 



Cf 



(3.14) 



All the results in this paper correspond to a fiducial 
cosmology in accord with the third year WMAP data 
release, i.e. Spergel et al. [2l[: flat ACDM universe with 
ft m /i 2 =0.128, n b ft 2 =0.02, /i=0.73, n s =0.96, r=0.089 and 
cr 8 =0.76. 

Figure [1] shows the variance of the residual 1Z(1), scaled 
by the CMB projected potential power spectrum, for the 
three examples that we consider. The SNAP (LSST) es- 
timator removes up to 90% (80%) of the CMB projected 
power for the lowest multipoles (1< 100) and only up to 
20% (10%) for K1000. This happens largely for two rea- 
sons. First and most important, LSST and SNAP do 
not have source galaxies at redshifts higher than «4 and 
«6, respectively. Therefore the potential fluctuations be- 
tween the upper- limit redshifts 4(6) and 1100 cannot be 
reconstructed; since the horizon is much smaller at these 
redshifts than today, the small scale power suffers greatly 
from this effect. Second, the shot noise starts to domi- 
nate the lensing signal at 1~1000, so the reconstruction 
is bound to fail at high multipoles, independent of the 
cosmological information provided by the source bins. It 
is easier to separate the impact of shot noise on the esti- 
mator from that of the redshift depth of the WL survey 



if we look at our third fiducial survey, the box distri- 
bution. The green, solid line shows the performance of 
the optimal estimator using the box distribution with 
Zmax =5 and n =100 galaxies/arcmin 2 . The total shot 
noise is the same as for SNAP, but the reconstruction is 
improved because there are more high-redshift sources. 
If we take z max =20 and keep the same total angular 
concentration of galaxies, then the estimator improves 
by a factor of «3 at low multipoles. We expect that by 
taking z rnax — > z CMB and n — > oo, we can recover entirely 
the projected potential to the CMB. The red solid curve 
in Fig. [TJ correponding to a redshift depth of 100 for the 
box distribution, seems to confirm this expectation. The 
estimator's efficiency increases relatively smoothly with 
the value of n : the difference between having 100 and 
250 galaxies/arcmin 2 leads to at most a factor of 2 in 
the estimator's efficiency for z max =100. This difference 
becomes smaller with lower z max . We conclude that for 
the purpose of the CMB projected potential reconstruc- 
tion, going to high redshift in a WL survey is much more 
important than trying to resolve a high number density 
of source galaxies. In section ijlVI we give a more quan- 
titative flavor to this conclusion and discuss what values 
of n and z max we need in a WL survey to significantly 
delens a _B-mode map. 



IV. DETECTING THE TENSOR-TO-SCALAR 
RATIO 



A. The Null-Hypothesis Test 



We would like to answer the following question: given 
a CMB survey, how much can we improve the detectabil- 
ity of the tensor-to-scalar ratio r if we use our projected- 
potential estimator to delens the CMB? There are sev- 
eral factors that influence the answer to this question: 
the sensitivity and angular resolution of the CMB exper- 
iment that measures the B modes, the noise level and 
redshift depth of the WL survey that provides the infor- 
mation for our potential estimator and, finally, the true 
value of r itself. The statistical significance of detection 
for some r > is estimated from the \ 2 statistic that 
would be obtained under the null hypothesis r — 0. The 
underlying approximation is that the delensed B mode 
obeys Gaussian statistics, which would be true for per- 
fect delensing and no foregrounds. We choose the sce- 
nario r = to be our fiducial model. We compare other 
models, with r ^ to the fiducial one and we determine 
the minimum r statistically distinguishable from r = 0. 
The delensed B del {l) will be constructed by subtracting 
from the lensed B(l) an estimate of the two lensing terms 
in equation (|2.ip . Accounting for the fact that we do not 
have access to the real lensing potential and E mode, but 
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FIG. 1: The residual variance of the CMB lensing poten- 
tial relative to its original value. LSST is the red dotted 
line, SNAP is the blue dashed line and the solid lines 
are for the box distribution. For the latter case, n =100 
galaxies/arcmin 2 and the thickness of the lines increases 
with the redshift depth of the survey. 




FIG. 2: The power spectra of the tensor B mode, the 
WL B-mode contamination and its residual after delens- 
ing. The black dotted and thin solid lines are the tensor 
B mode power spectra for r = 10~ 4 and r = 10~ 2 , respec- 
tively. The upper solid black line is the full WL B mode 
power and the rest of the lines are the residual contam- 
ination power spectra after delensing with our optimal 
estimator. The latter correspond to u>~ 1//2 = ar- 
cmin and #fwhm=10 arcmin. The red, long-dashed curve 
corresponds to # F whm=30 arcmin and z ma x =100. 



their estimators, the delensed field is 
d 2 V 



B del {l) = B{1) 



2tt 

1 f <Pl'<Pl" 

2 2vr 2tt 



f(l',l) [HL)E(l') - i>{L)E{V) 
g(l',l",l)[TP(l")r(l"-L)E(l>) 



-ip(l")ip*(l"— L)E(l') ] + B (I). 



(4.1) 



We include the measurement noise B N , whose variance is 
given by equation ()3 . 1 3[) . The variance of the measured, 
delensed B mode follows from the above: 



(jB, del ^^jtensor 



cf 1 



WL, del 



We can express the residual power left from delensing, 
£<wL,dei ^ j n -j- erms Q £ tfie variance of the residuals of E 

and ip, and the power spectra of the real lensing potential 
and E mode: 



C, 



wl, del 



The last integral in the above equation, containing terms 
of the second order in the potential power, yields a much 
smaller contribution to (j™ L - del compared to the first 
term. For computational simplicity, we choose therefore 
to neglect this last term; all the results following have 
been obtained in this approximation. Let us write the 
Ax 2 for the measured and delensed B-mode power spec- 
trum: 



A X 2 {r)=Y[c^ dcl {r)-C^ dd {r = 0) 



Var 



s-iB, del / \ 

Ci (0 



(4.4) 

(4.2) The variance of the delensed B-mode power spectrum is: 



Var(Cf dei (r)) 



(21 + l)f sky 



-Cf 



C, N + C, 



wl, del 



-Var^Var^] + i|^^| 5 (r,r,0| 2 x 
{Cf„ C r _ L \ YwTlf + ( - Varft^ ) [ C% Varftj;,, + 
+C\i"- L \ V&rTlf,, - Var nf„ VaxTZt,_ L , ] }. (4.3) 



(4.5) 

Here f s k y denotes the fraction of sky covered by the CMB 
mission. To be able to use our lensing potential estima- 
tor, the CMB experiment must survey the same area of 
the sky as the WL experiment. Using the above equa- 
tion in the case of the fiducial model r = 0, we rewrite 
Eq. JOJ 



as: 



Ax 2 (r)=£ 



21- 



J sky 



j, ^jtensor 
£,WL,del Qjy 



(4.6) 
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In order to find the smallest r discernible from the fiducial 
value of 0, we need to solve the following equation: 

A X 2 (r) = a 2 , (4.7) 

where a sets the confidence level. All our results are com- 
puted for a— 2, which corresponds to a 95% confidence 
level for a single free parameter r. 

B. Results 

Figure shows the lensing B mode, given by Eq. 12.61 
and the tensor B mode for two values of r = 10~ 2 , 10~ 4 . 
The latter was computed using the publicly available 
code CAMB, i.e. [3(|, Lewis et al. [Hj]. Also shown is the 
residual lensing power, given by Eq. 14.31 when delensing 
using LSST and box estimators has been applied. With 
one exception, the plots correspond to a beam $fwhm = 10 
arcmin; for all of them, the noise level is w^ 1 / 2 = 1/iK 
arcmin, i.e the expected CMBPol noise level. 

We can see that LSST delensing can reduce the lens- 
ing B-mode power by a factor of wl.3 («2 for SNAP, 
although not shown in the figure), while a box with 
Zmax =20 and n =100 galaxies/arcmin 2 leaves a residual 
power «8 times lower than the original contamination. 
Although such high redshifts are not attainable by a WL 
galaxy survey, we would like to illustrate the potential 
benefits of a 21-cm emission survey. There is further 
improvement as we go to higher z max : for z max =100 
and the same value of n , the residual lensing power is 
«37 times lower than the full lensing .B-mode power. In- 
creasing the beam width has a dramatic effect: the long- 
dashed red line is for the same box with z max =100, but a 
beam with # FWHM =30 arcmin; compared to the 10 arcmin 
beam, the delensing efficiency has dropped by a factor of 
6. This is understood as follows: the lensed -B-mode 
at I ~10-100 arises in part from the beating of -E-mode 
power at I ~ 1000 with deflection power at similar I. The 
delensing must know both these E and deflection fields at 
I « 1000 in order to succeed, hence the CMB experiment 
must resolve these modes. 

Naturally, we would like to probe now how much de- 
lensing with the above-mentioned WL survey models im- 
pacts the detection of tensor modes. In figure [3] we plot 
the minimum tensor-to-scalar ratio distinguishable from 
at 95% confidence level as a function of the CMB exper- 
imental noise for a fixed beam (9 FWHM =10 arcmin. r m i n 
is obtained by solving Eq. 14.71 for some of the delens- 
ing scenarios discussed earlier and the noise is in the 
form of weight per solid angle, w~ x / 2 . The upper thin, 
solid, black curve corresponds to the case when no de- 
lensing is done. We also show (turquoise and red arrows) 
the approximate noise levels w" 1 ^ 2 — 60 and 1 /iK ar- 
cmin for two future missions: the Planck satellite [3lJ 
and CMBPol. The noise level for the latter was es- 
timated by assuming 1000 detectors with a sensitivity 
of 50 /iK^/sec, sky coverage of 20000 deg 2 and 5 years 
of operation. Delensing is relevant for the detection of 



tensor B modes only if the instrumental noise is low 
enough, to -1 ' 2 < 10/iK arcmin. And this is an opti- 
mistic value, since our analysis does not take foregrounds 
into account. While CMBPol is meant to reach such low 
noise levels, delensing will not, most probably, be a con- 
cern for Planck. This is in accord with the findings of 
the Planck Science Team. For low noise levels, when de- 
lensing does play a part, figure [3] reinforces our earlier 
conclusions: the difference between the case of no de- 
lensing and that of a box distribution going to z max =20 
translates approximately into a factor of 7 difference in 
the corresponding r m i n for the lowest detector noise that 
we consider, w~ x / 2 = 0.1/iK arcmin. z max =100 gains 
another order of magnitude in the detection threshold of 
the tensor B modes, for the same weight per solid angle. 
This is quite idealized, however: for the technically ambi- 
tious CMBPol, at w~ 1/2 = 1/xK arcmin, the extension to 
Zmax =100 gains only a factor of 3 over z max =20. The 
21-cm signal does not exist significantly above z k 100 
since the hydrogen spin temperature is expected to match 
the CMB temperature at earlier epochs, for instance see 
Barkana and Loeb [23| and Zahn and Zaldarriaga 24]. 

We note that r min corresponding to LSST delensing is 
lower by a factor of «3 than the SNAP r m i n , despite the 
better performance of the SNAP lensing potential esti- 
mator. This is due to LSST's superior sky coverage and 
the presumed larger CMB sky coverage. We emphasize 
that the improvement in r m i n for LSST and SNAP de- 
lensing relative to the case of no delensing (correspond- 
ing to their respective survey areas of 20000 deg 2 and 
1000 deg 2 ) is of only w20% and 50%, respectively. We 
expected this result: the potential modes with angular 
scales of I wlOOO mix with E power and generate the 
low-Z B mode polarization power; but such high-/ lensing 
potential modes are poorly reconstructed by both LSST 
and SNAP. 

Figure 01 illustrates the importance of the beam choice 
for our delensing method. For the CMBPol noise level 
w^ 1 / 2 = 1/tK arcmin, we plot r m i n as a function of beam 
size, for the same box models mentioned in Fig. [31 If 
fl FWEM <10 arcmin, r min stays constant for all models 
considered. But once 6* FWHM >10 arcmin, the delens- 
ing efficiency drops steeply and for # FW hm = 2°, there is 
nothing to be gained from delensing. 

We now explore the importance of the WL free pa- 
rameters z max and h for the tensor B modes detection. 
Figure [5] is a contour plot of r m i n as a function of red- 
shift depth and angular concentration of source galaxies. 
All contours are plotted using the box distribution for a 
CMBPol-level of detector noise, ^ FWHM — 10 arcmin and 
a 20000 deg 2 sky coverage. This plot confirms our find- 
ings from section SpTTCj Once a certain source density is 
reached (n «50 galaxies/arcmin 2 at z max <10), there is 
no significant improvement in r m i n if we continue to in- 
crease n . The redshift depth is the dominant parameter 
of the WL survey: even with a concentration as low as 10 
galaxies/arcmin 2 , if we reach sources at z=15, r m i n has 
a factor of 2 improvement over the case of no delensing. 
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r. 




FIG. 3: The minimum tensor-to-scalar r distinguishable 
from at 95% confidence level as a function of CMB de- 
tector noise. All curves are for n = 100 galaxies/arcmin 2 , 
with the exception of the LSST line, for which n =30 
galaxies/arcmin 2 . The upper blue dashed curve is for 
SNAP delensing (.4=1000 deg 2 ). The rest of the lines 
correspond to a 20000 deg 2 sky coverage and match the 
delensing cases presented in Fig. [2] We have also indi- 
cated the approximate polarization noise levels for the 
Planck and CMBPol missions. 
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FIG. 4: Dependence of r m i n on the beam size. The curves 
are for a box distribution with z max =5, 10, 20 and 100, 
n =100 galaxies/arcmin 2 and the noise level of CMBPol. 



Only if we reach z max >100 can we significantly improve 
Train by increasing n above 50 galaxies/arcmin 2 . 



C. Related Work in the Literature 

Knox and Song [9( investigate the detectability of in- 
flationary gravitational waves using CMB observables to 
delens. They apply the quadratic estimator of Hu and 
Okamoto Q to a CMB polarization experiment and thus 
reconstruct the lensing potential. They find that for a 
noise-free CMB map, the lensing .B-mode power can be 
reduced by an order of magnitude, which corresponds to 
a similar reduction in r m i n . 

Kesden et al. [8j apply the same quadratic estimator to 
temperature maps and find similar results as Knox and 
Song 0. 

Sigurdson and Cooray [l4j reconstruct the lensing po- 
tential using the above-mentioned quadratic estimator 
applied to 21-cm temperature fluctuations. They account 
for two distinct sources of noise that affect their estima- 
tor: the reconstruction noise, and the unreconstructed 
lensing due to mass between the 21-cm source redshift 
and the CMB. In the limit where the reconstruction noise 
dominates, they find that a 21-cm survey, combined with 
a very sensitive CMB survey, could reduce the detection 



threshold of the tensor i?-mode power by an order of 
magnitude, compared to the no delensing case. The au- 
thors also consider the limit where the partial delensing 
bias dominates. In this limit, the reduction in the lensing 
contamination power is the same as what we find if we 
take n — > oo. 

Seljak and Hirata [25j show that the reconstruction of 
the lensing potential using the maximum likelihood esti- 
mator applied to polarization maps can reduce the power 
of the lensing B mode by at least a factor of 40. They 
conclude that delensing is not a fundamental limit to 
the detection of inflationary gravitational waves, if high- 
resolution CMB maps cleaned of foreground emission are 
available. 



V. CONCLUSIONS 

In this paper we have examined the possibility of de- 
lensing i?-mode polarization maps using galaxy WL sur- 
veys. We have proposed a weighted combination of pro- 
jected potential estimators to different source redshift 
bins which optimally reconstructs the projected poten- 
tial seen by the CMB. We have used three fiducial sur- 
veys to exemplify our estimator: LSST, SNAP and a 
generic survey relevant mostly to future 21-cm studies. 
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FIG. 5: r min as a function of the WL survey parameters, 
Zmax and n. The contours are for a box distribution, the 
noise level of CMBPol and # FWH m = 10 arcmin. Their increas- 
ing thickness depicts the change in the detectable tensor B 
mode amplitude as we transit from a weak delensing case 
{zmax =1, n =250 galaxies/arcmin 2 ) to very good delensing 
{zmax =100, n =250 galaxies/arcmin 2 ). 

These examples have different source redshift distribu- 
tion, source density n , and sky coverage / s ky, and have 
enabled us to test the effect of each of these factors on 
the lcnsing potential estimator and also on its ability to 
reduce the WL contamination of B mode maps. Using 
a Ax 2 test for the delensed B mode field, we have de- 
termined the minimum value of the tensor-to-scalar ratio 
statistically distinguishable from in optimally delensed 
data. Throughout this paper we have ignored the po- 
larized foreground contamination and we have also made 
the assumption of Gaussianity for the delensed B mode 
map. While this assumption may be inaccurate, espe- 
cially for the cases when the efficiency of delensing is low, 
we expect the qualitative results of our work to hold. 

The lensing potential estimator is sensitive mostly to 
the redshift depth of the WL survey and reconstructs 
best the large angular scale multipoles. The performance 
of the estimator improves continually as higher redshift 
source galaxies are used. In the limit z rnax —> z OMB and 
n — ► oo, the reconstruction is perfect. An experiment like 
SNAP recovers «90% of the CMB projected potential 
power for I < 100 and «20% for I < 1000. The perfor- 
mance of LSST is a little worse, because it detects sources 
at an average redshift of {z) = l, compared to (z)=1.5 of 
SNAP. A box distribution (constant redshift distribution 
of source galaxies), provides better CMB potential esti- 
mator: if we go as far as z max =20, there is an order of 
magnitude improvement over SNAP, for the same angu- 



lar concentration of galaxies. 

However, in this last case, the reduction in r m i n , com- 
pared to the case where no delensing is done, is only of 
a factor of «7 for u>~ 1//2 = 0.1/iK arcmin and # FWH m = 10 
arcmin. This happens because low-/ lensing B-modes are 
generated by beating of gravitational potential modes 
and B-modes on scales coresponding to ZwlOOO, where 
the potential estimator is less faithful. In the case of 
LSST and SNAP, the reduction in r m i n relative to the 
no delensing case is only by «15% and «50% for the 
same CMB noise level and beam. 

We stress that delensing with foreground weak lensing 
offers significant gains in r m i n (factors of 5 or more) only 
when three conditions are met: 

1. the noise in the CMB map is sufficiently low, ~ 
1 [iK arcmin; 

2. the beam size of the CMB map is < 20 arcmin; 

3. the lensing source distribution extends to z ~ 20 
or greater. 

Delensing is not relevant to tensor mode detection if 
w- 1 / 2 > 10/iK arcmin or if FWHU >2°. Also, there is no 
advantage in lowering the beam size beyond # FWH m = 10 
arcmin: for our delensing method, a 1 arcmin beam will 
do just as well as a 10 arcmin beam. This is perhaps 
the most important feature of our delensing technique, as 
both the quadratic estimator of Hu Q and the maximum 
likelihood estimator proposed by Hirata and Seljak [ll| 
require beam sizes of 2-3 arcmin to yield their best per- 
formance. For the CMBPol mission, (which might reach 
a noise of 1 /iK arcmin) delensing with a box of z max =20 
results in r m i n « 2 x 10~ 5 . If no delensing is applied, then 
Train ~ 8 x 10 -5 , a factor of 4 worse. CMBPol detector 
noise keeps r m i n > 7 x 10~ 6 even with perfect delens- 
ing, so there is less incentive to acquire delensing source 
planes above z max =20. Also for the CMBPol noise level, 
we have investigated the impact of n and have concluded 
that as long as we go to high redshift (z max > 10), even 
with a low density of source galaxies, we can still improve 
Tmin by a factor of a few. 
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